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Abstract 



J^ ■ We calculate the relic density of dark matter in the MSSM with CP violation. We 

ly-^ , analyse various scenarios of neutralino annihilation: the cases of a bino, bino-wino 

and bino-Higgsino LSP, annihilation through Higgs, as well as sfermion coannihi- 
lation scenarios. Large phase effects are found, on the one hand due to shifts in 
\^ , the masses, on the other hand due to modifications of the couplings. Taking spe- 

^^ I cial care to disentangle the effects in masses and couplings, we demonstrate that 

^S^ • the presence of CP phases can have a significant influence on the neutralino relic 

'j^. abundance. Typical variations in Qh"^ solely from modifications in the couplings are 

Qh! C'(10%-100%), but can reach an order of magnitude in some cases. 

^^ 
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1 Introduction 

With the conclusive evidence for a significant component of cold dark matter (CDM) in 
the Universe, there is considerable interest, both at the theoretical and experimental level, 
to identify this CDM and analyse its properties; see [1] for a recent review. In particular, 
if the CDM consists of a new weakly interacting massive particle (WIMP), as predicted 
in generic new physics models with a discrete symmetry that ensures the stability of 
the lightest particle, the next generation of colliders has good prospects to discover it. 
Being electrically neutral and stable, the WIMP escapes the detector as missing energy 
and momentum. The preferred discovery channels therefore rely on the production of 
other new particles present in the theory and their decays into the CDM candidate. By 
measuring the properties and decay kinematics of these new particles, one should then be 
able to determine the properties of the WIMP. If the measurements are precise enough, 
this allows to predict the annihilation cross sections and hence the thermal relic density 
of the CDM candidate, thus checking the consistency between a particular model of new 
physics and cosmology. 

Indeed, particle physics models trying to explain the dark matter are constrained 
by recent precision cosmological measurements. These are most notably the data from 
WMAP [2,3] and SDSS [4], which imply a (dominantly) cold dark matter density of 

0.0945 < ncDuh^ < 0.1287 (1) 

at 2cr. In the following we will refer to Eq. (1) as the WMAP range. 

The relic density of dark matter has been discussed extensively in the framework of 
the most popular model for new physics, low-scale supersymmetry (SUSY) with R-parity 
conservation. Especially if the lightest supersymmetric particle (LSP) is the lightest neu- 
tralino, this provides a good cold dark matter candidate [5,6]. For the standard picture 
of thermal freeze-out [7], assuming no additional non-thermal production mechanism, 
comprehensive public codes that compute the relic density of the neutralino LSP are 
available today: micrOMEGAs [8,9], DarkSUSY [10], and IsaRED [11]. Many analyses of 
neutralino dark matter were performed in the minimal supersymmetric standard model 
(MSSM), whose parameters are defined at the weak scale, see e.g. [12-14], in the con- 
strained MSSM (CMSSM) and mSUGRA models [15-25], and in other models with the 
parameters defined at the unification scale [26-34]. 



These studies showed that in the MSSM there are only a few mechanisms that provide 
the correct amount of neutrahno annihilation, consistent with the WMAP range Eq. (1): 
annihilation of a bino LSP into fermion pairs through t-channel sfermion exchange in case 
of very light sparticles; annihilation of a mixed bino-Higgsino or bino-wino LSP into gauge 
boson pairs through t-channel chargino and neutralino exchange, and into top-quark pairs 
through s-channel Z exchange; and finally annihilation near a Higgs resonance (the so- 
called Higgs funnel). Furthermore, coannihilation processes with sparticles that are close 
in mass with the LSP may bring r2/i^ in the desired range. In particular, coannihilation 
with light sfermions can help to reduce the relic density of a bino-like LSP. Here note that 
coannihilation generically occurs when there is a small mass gap between the LSP and 
the next-to-lightest SUSY particle (NLSP). In scenarios with a Higgsino or wino LSP, 
one has in fact a mass-degenerate triplet of Higgsinos or winos, and coannihilations are 
so efficient that ^h? turns out much too small, unless the LSP has a mass of order TeV. 
In the case that r2/i^ is too low one would need, for instance, a significant contribution 
from non-thermal production. 

Most of these previous analyses assumed that CP is conserved, although CP- violating 
phases are generic in the MSSM. Here note that, given the Higgs mass bound ruh > 
114 GeV from LEP [35], the CKM description of CP violation (CPV) in the SM is not 
sufficient to generate the correct baryon asymmetry in the Universe. One possible solution 
to this problem is leptogenesis. Another solution is electroweak baryogenesis through new 
sources of CP violation beyond the SM. Electroweak baryogenesis in the MSSM with CP 
violation (CPV-MSSM) has been studied in [36-38]; see also references therein. For a 
strongly first-order phase transition, it relies on the existence of a light stop [39-42] . 

The parameters that can have CP phases are the gaugino and Higgsino mass param- 
eters and the trilinear sfermion- Higgs couplings. Although constrained by electric dipole 
moments, nonzero phases can significantly infiuence the phenomenology of SUSY parti- 
cles. They can also have a strong impact on the Higgs sector, inducing scalar-pseudoscalar 
mixing through loop effects [43,44]. This can in turn have a potentially dramatic effect on 
the relic density prediction in the Higgs funnel region: neutralino annnihilation through 
s-channel scalar exchange is p-wave suppressed; at small velocities it is dominated by 
pseudoscalar exchange. In the presence of phases, both heavy Higgs bosons can, for 
instance, acquire a pseudoscalar component and hence significantly contribute to neu- 
tralino annihilation even at small velocities. Likewise, when only one of the resonances 
is accessible to the neutralino annihilation, large effects can be expected by changing 
the scalar /pseudoscalar content of this resonance. Last but not least, the couplings of 
the LSP to other sparticles depend on the phases, and so will all the annihilation and 
coannihilation cross sections, even though this is not a CP-violating (CP-odd) effect. 

Consequences of complex parameters for the relic density of neutralinos have so far 
been considered in the literature only for specific cases. The effect of a phase in the squark 
left -right mixing was shown to enhance the cross section for neutralino annihilation into 
fermion pairs [45], but the important sfermion coannihilation effects were not included. 
In [46], the influence of phases was discussed for various cases gaugino- Higgsino mixing, 
taking into account CPV in the Higgs sector. However, coannihilations were again dis- 
regarded. Reference [47] discussed phase effects in neutralino annihilation through Higgs 
exchange in the CMSSM at large tan/3. Likewise, Ref. [48] examined a very speciflc case, 
namely the dependence of Qh"^ on phases that arise through supersymmetric loop correc- 
tions to the b quark mass in SUGRA models with Yukawa uniflcation. For the general 



MSSM, the effect of a phase in the trihnear couphng of the top squark on the neutrahno 
annihilation near a Higgs resonance were discussed in [49] and more recently in [50] , how- 
ever without analysing the phase dependence of Qh^. Finally, the case of coannihilation 
with very light stops was discussed in [36,38] in the framework of electroweak baryogene- 
sis. Also these latter papers considered only few fixed values of 0^ and did not study the 
phase dependence of Qh"^ in detail. It is, moreover, important to note that all the above 
mentioned publications studied the impact of phases solely in terms of SUSY-breaking 
parameters, without discussing that part of the phase effects are due to changes in the 
masses of the involved particles. 

A complete and coherent analysis of the relic density of neutralinos in the CPV- 
MSSM is therefore still missing, and this is the gap we intend to fill with this paper. 
We perform a comprehensive analysis of the impact of CP phases on the relic density of 
neutralino dark matter, taking into account consistently phases in ALL annihilation and 
coannihilation channels. To this end, we realized the implementation of the CPV-MSSM 
within micrOMEGAs 2.0 [51,52]. For the computation of masses, mixings and effective 
couplings in the Higgs sector, we rely on CPsuperH [53]. We further take into account 
collider constraints from sparticle and Higgs searches, as well as the constraint arising 
from the electric dipole moment of the electron [54] . 

Besides realising the complete calculation of the relic density of neutralinos in the CPV- 
MSSM, we analyse in detail the above-mentioned scenarios of neutralino annihilation and 
coanihilation, for which the LSP is a 'good' CDM candidate. We find indeed a large 
impact from phases due to modifications in the sparticle couplings (through changes in 
the mixing matrices) but also due to changes in the physical masses. Some of the largest 
effects are in fact due to kinematics. This should be expected as the relic density is 
often very sensitive to masses, in particular to the exact mass difference between the LSP 
and NLSP in coannihilation processes, or in the case of annihilation near a s-channel 
Higgs resonance, to the difference between twice the LSP mass and the mass of the 
Higgs. In these scenarios, setting apart the purely kinematic effects hence somewhat 
tames the huge effects due to CP phases found in some of the previous studies listed 
above. On the other hand, we also find cases where the phase dependences of masses and 
couplings work against each other, and taking out the kinematic effects actually enhances 
the phase dependence of the number density. There are even cases in which, for fixed 
MSSM parameters, Qh"^ decreases with increasing phase, but once the masses are kept 
fixed, Qh"^ actually grows. Since what are relevant to experiments are rather the physical 
observables (masses, branching ratios, etc.) than the underlying parameters, we take 
special care in our analysis to disentangle effects arising from changes in the couplings 
from purely kinematic effects. 

After a description of the model in section 2 and its implementation into micrOMEGAs 2 . 
in section 3, we present in section 4 our results for the typical scenarios of neutralino 
annihilation: the mixed bino-Higgsino case with annihilation into W pairs, rapid an- 
nihilation near a Higgs resonance, the bino-like LSP with light sfermions (t-channel and 
coannihilation regions), and finally the mixed bino-wino scenario with neutralino-chargino 
coannihilation. In section 5, we give a summary and conclusions. For completeness, the 
Lagrangian for sparticle interactions in the CPV-MSSM is given in the Appendix. 



2 The MSSM with CP phases 

We consider the general MSSM with parameters defined at the weak scale. In this frame- 
work, the gaugino and Higgsino mass parameters and the trilinear couplings can have 
complex phases, Mi = \Mi\e'^'^', /i = |yu|e*'^'" and Aj = \Af\e^'f'f . Not all of these phases are 
however relevant to our analysis. In particular, the phase of M2 can always be rotated 
away, while the phase of M3 is mostly relevant for the coloured sector, so we neglect it in 
this study. Among the trilinear couplings. At has the largest effect on the Higgs sector, 
and can thus potentially play an important role for the relic density. The phase of the 
selectron coupling A^, as the phase of all light sfermions, is irrelevant for the relic density; 
it has, however, to be taken into account since it contributes to electric dipole moments 
(EDMs). 

In the following, we explain our notation and conventions. We hereby basically follow 
the notation in CPsuperH and use the SUSY Les Houches Accord (SLHA) [55]. We also 
discuss some of the most relevant couplings. For the complete interaction Lagrangian, see 
the Appendix and Ref. [53]. 

Neutralinos: The neutralino mass matrix in the bino-wino-Higgsino basis 

^o = HA', -^3,^^° ,7^0 )i 



IS 



M 



N 



I Ml -mzswcp mzSwSfB \ 

M2 mzCwC/3 -mzCwSfs 

—mzSwCfB nizcwcp — /x 

\ mzswsp -nizCwSfs -/i / 



(2) 



with sw = sm9]v, cw = cos9w, s^ = sin/5, cp = cos/3 and tan/5 = 1^2/^1 (^1,2 being the 
vacuum expectation values of the two Higgs fields -^1,2)- This matrix is diagonalized by 
a unitary mixing matrix A^, 

N*A4j\iN'' = diag(m^o, m^o, m^o, m^o) , (3) 

where m^o, i = 1,...,4, are the (non- negative) masses of the physical neutralino states 
with m^o < .... < m^o. The lightest neutralino is then decomposed as 

X? = NuB + Ni.W + NrsHi + A^i4i?2 (4) 

with the bino (Jb), wino (fw) cind Higgsino (fn) fractions defined as 

fB = \Nu\\ fw = \Ni2\\ fH = \N,3\^ + \N,4\ (5) 

The LSP will hence be mostly bino, wino, or Higgsino according to the smallest mass 
parameter in Eq. (2), Mi, M2, or /i, respectively. 

Charginos: The chargino mass matrix in SLHA notation, ^ 

Mc=( ^ ^' v/2m^sin/3\ 

\ ^/zniw cos/5 jj, J 



^Here note that in CPsuperH, Mc -^ M 



T 

c- 



is diagonalized by two unitary matrices U and V, 

U*McV^ = di^g{m^±,m^±), (7) 

with the eigenvalues again being mass-ordered. The chargino sector will be relevant mainly 
because t-channel chargino exchange dominates the XiXi ~^ W^W~ cross-section, which 
is often the main annihilation process when the LSP has a sizeable Higgsino component. 
The relevant Lagrangian writes 

^x?x+VF- = 9W;^^^^ {Of.PL + Of.Pn) xt + h-c. (8) 

with 

o^ = N,2V^ - -^N,,v:, , Of, = N;,Ua + -^n;,u,2 . (9) 

These expressions are the same as for the real MSSM and show that the XiXj W^^ cou- 
pling vanishes in case of a pure bino LSP. Note that this coupling also enters the chargino 
coannihilation processes into gauge boson pairs which are important in some specific sce- 
narios. Nonzero phases will affect directly the masses of the charginos and neutralinos as 
well as the couplings through shifts in the mixing matrices. 

Sfermions: Ignoring intergenerational mixing, the sfermion mass matrices are 

tTTl h* fl* 

. ^" ^2^ ] (10) 



with 



^2 _ „^2 , ^2 ^^^o/QfTf _„.ci^2^...^ I ^2 



j^ - Mf^j^^ + m'z cos2(3 m^-ef sin' ew)+mj, (11) 

2 = M^ 



m] = M'^jy D R| + *^/ "^1 cos 2(3 sin^ 6w + rrij, (12) 



and 



at = {Av2-fi*v,)/V2, (13) 

ab,r = {Ab,rVi-fi*V2)/V2. (14) 

for f = t, b, f. Here Mrg i^, ^w dr\ ^^^ ^^^ ^^^ ^"^^ right sfermion mass parameters; 
mj, 6/ and I^ are the mass, electric charge and the third component of the weak isospin 
of the partner fermion, respectively, and hf denote the Yukawa couplings, ht = ^^0 and 

hhT = — ^-^ at tree level. Here, w^ = w? + v^- The mass matrix A^^- is diagonalized by a 
unitary matrix R^ such that 

R^'^ M)R^ = diag(mj-^, mp (15) 

with m? < m? and (/^, f^)^ = R^{f-j^, j^)'^- The relevant sfermion interactions are 
given in the Appendix. 



Neutral Higgs bosons: In the CP-conserving MSSM, the Higgs sector consists of two 
CP-even states h^, H^ and one CP-odd state A^. Allowing for CP- violating phases in- 
duces a mixing between these three states through loop corrections. The resulting mass 
eigenstates /ii, /i2, h^ (with rrih^ < rrih^ < nih^ by convention) are no longer eigenstates of 
CP. The Higgs mixing matrix is defined by ^ 

(01,02,a)J = ifai(/il,/i2,/i3)f- (16) 

Because of the mixing between the neutral scalar and pseudoscalar states, it is preferable 
to use the charged Higgs mass, mfj+, as an independent parameter. In what follows we 
will mainly be concerned with the couplings of the lightest neutralino to Higgs bosons 
that govern the neutralino annihilation cross section via Higgs exchange. The Lagrangian 
for such interactions writes 

3 

^x?x^. = -f 5^ ^i(4x;x? + ^^s^CxSx?)*?^^ (1'^) 

with the scalar and pseudoscalar couplings corresponding to the real and imaginary part 
of the same expression: 



4x?x; = Re[{N*,-twN*,){HuN*,-H,,N*,-tHs.{spN*,-Cf,N*,))], (18) 
<x?x? = -Im[{N;,-twN*,){HuN;,-H,,N*,-tHs.{sfsN*,-Cf,N*,))]. (19) 



For real parameters, one recovers the {h^, H^, A'^) couplings by taking for the pseudoscalar, 
A^, H^^ = 1 and all other Hij = 0, while for the scalar Higgs couplings, the only nonzero 
elements are Hu = —H22 = — sin a and -^21 = ii\2 = cos a. The LSP couplings to neutral 
Higgs bosons will clearly be affected both by phases in the neutralino sector which modify 
the neutralino mixing, as well as by phases which induce scalar-pseudoscalar mixing in the 
Higgs sector. Indeed in the MSSM the Higgs CP mixing is induced by loops involving top 
squarks and is proportional to Im{Atfi) / [m'^ ~^^'i ) t^^]. Thus a large mixing is expected 
when Im{At^) is comparable to the stop masses-squared. Note that the masses of the 
physical Higgs bosons also depend on the phases of At and /i. In particular, a large mass 
splitting between the heavy states is found for large values of jjiAt. For a more detailed 
discussion of Higgs masses and couplings in the CPV-MSSM, we refer to 



3 Relic density computation and implementation of 
the CPV-MSSM into micrOMEGAs 

The computation of the relic density of dark matter is by now standard [57,58] . It relies on 
solving the evolution equation for the abundance, defined as the number density divided 
by the entropy density. 



^ = V ^^^^ <^^') (^(^)' - ^^^(^)') = ^(^) (^(^)' - ^^«(^)') ' (20) 



^In CPsupcrH, the Higgs mixing matrix is denoted Oai 



where g^ is an effective number of degrees of freedom [57], Mp is tlie Planck mass and 
Y(,q{T) tlie tliermal equilibrium abundance, and {av) the relativistic thermally averaged 
annihilation cross-section of superparticles summed over all channels. It is given by 

, . »J (m,+m,)2 k,l 

(av) = 7j , (21) 

2T{Z9^m^K2{m,/T)Y 

i 

with Qi the number of degrees of freedom, aij-ki the total cross section for annihilation 
of a pair of supersymmetric particles with masses rrii, rrij into some Standard Model 
particles {k,l), and pij (y/s) the momentum (total energy) of the incoming particles in 
their center-of-mass frame. Integrating Eq. (20) from T = oo to T = Tq gives the present 
day abundance Y(Tq) needed in the evaluation of the relic density, 

n^sph' = ^^^^MLSPr(ro) = 2.742 X W'^Y{n) (22) 

Per <jeV 

where s{Tq) is the entropy density at present time and h the normalized Hubble constant. 
The present-day energy density is then simply expressed as plsp = 10.54r2/i^(GeV/m^). 
To compute the relic density, micrOMEGAs solves the equation for the abundance, 
Eq. (20), numerically without any approximation. In addition, micrOMEGAs also estimates 
the relative contribution of each individual annihilation or coannihilation channel to the 
relic density. For this specific purpose, the freeze-out approximation is used: below the 
freeze-out temperature T/, Fgg -C Y, and Eq. (20) can be integrated 



A{T)dT. (23) 

To 



F(To) Y{Tf) 



It turns out that Mi^sp/Tf ~ 25 (for more details see [8,57]). In this case, the first 
term in Eq. (23) is suppressed, and one can approximate 1/Q by a sum over the different 
annihilation channels: 



t,j,k,l 



VlK^ .^, uJij-kih"^ 



where 



P-r / V 45 ' 2r(E9.™?/f2(m,/r)) = 

Note that this method gives to a good approximation the contribution of individual 
channels, although the value of Vth? is slightly overestimated as compared to the exact 
value obtained by solving directly the evolution equation for the abundance as described 
above. 

To perform this computation for the CPV-MSSM, we are using micrOMEGAs 2.0 [51], 
an extension of micrOMEGAsl .3 [8,9] that allows the computation of the relic density 



of dark matter within generic particle physics models that feature a stable weakly in- 
teracting particle. Within this framework we have implemented the CPV-MSSM as fol- 
lows. Using LanHEP [59], a new CPV-MSSM model file with complex parameters was 
rebuilt in the CalcHEP [60] notation, thus specifying all relevant Feynman rules. For the 
Higgs sector, an effective potential was written in order to include in a consistent way 
higher-order effects [44]. Masses, mixing matrices and parameters of the effective po- 
tential are read directly from CPsuperH [53], together with masses and mixing matrices 
of neutralinos, charginos and third generation sfermions. The masses of the first two 
generations of sfermions are computed at tree-level from the MSSM input parameters 
within micrOMEGAs. The cross sections for all annihilation and coannihilation processes 
are computed automatically with CalcHEP. The standard micrOMEGAs routines are used 
to calculate the effective annihilation cross section and the relic density of the LSP. This 
CPV-MSSM version of micrOMEGAs has first been presented in [52]. 

EDM constraint 

We have also implemented the constraints originating from the electric dipole moment 
of the electron (eEDM), 4 < 2.2 x IQ-^^ecm [54]. In the MSSM, for small to in- 
termediate values of tan/5, the dominant contribution to rfg originates from one- loop 
chargino/sneutrino and neutralino/selectron exchange diagrams [61,62]. At one loop, the 
expression for de depends on the complex parameters in the neutralino/chargino sector 
(Ml, /i) as well as on the trilinear coupling of the electron, Ag. Note that Ae is suppressed 
by a factor proportional to nif. and is completely negligible for cross sections calcula- 
tions. It is, however, relevant for the eEDM since rfe is itself proportional to the electron 
mass. The eEDM features a strong dependence on /i because the necessary helicity flip 
originates from the coupling of the electron to the Higgsino, which is proportional to /i. 
The one-loop contributions to the eEDM basically restrict 0^ to | sin0^| ~ (9(10~^) with 
sfermions at the TeV scale. The restrictions on the phases 0i and 0e are more modest. 
At one loop, the eEDM constraint can be most easily evaded by raising the masses of the 
sfermions of the first generation, say to 10 TeV. However, for such high sfermion masses, 
two-loop contributions [63,64] can also become important, especially if At,b and /i are also 
in the TeV range. Two-loop contributions are further enhanced for small values of the 
charged Higgs mass and for large tan/3 [65]. In the scenarios which we consider here, we 
will sometimes have to rely on a cancellation between one- and two-loop contributions to 
have an acceptable value for the eEDM. 

4 Results 

We now turn to the numerical analysis and present results for various scenarios for which 
the relic density is in agreeement with WMAP. We impose GUT relations for the gaugino 
masses. Mi : M2 : M3 ^1:2:6, unless mentioned otherwise. In order to satisfy 
the eEDM constraint, we assume in most cases that the sfermions of the first and second 
generation are heavy, m ? =10 TeV, allowing only the third generation to be at the TeV 
scale. Unless specific values are specified, we take Ms = Mq^ = Mjj^ = M^^ = Mi^ = 
Mj^ . For the trilinear couplings, we keep At as a free parameter, assuming \Af\ = 1 TeV 
for all others. This is justified as the Af with / 7^ t have only a very mild effect on the 
neutralino cross-sections. However, one has to keep track of A^ because the electron EDM 



depends on the phase in the slepton sector. For simphcity we consider a common phase 
(pi for all trilinear slepton couplings. In general, the parameters that will be allowed to 
vary hence are 

|Mi|, |/i|, tan/3, mH+, |A|, Mg, 0i, 0^, 0i, (f)i . (26) 

As mentioned above, the eEDM will constrain 0^ to ~ or 180°. From now on we drop 
the II for simplicity of notation, i.e. |Mi|e*'^i — »■ Mie*'^^, etc. So when specifying the value 
of a complex parameter, we implicetly mean its absolute value. 

The cross sections for the annihilation and coannihilation processes will depend on 
phases, and so will the thermally-averaged cross section {(tv), Eq. (21). Part of this is 
due to changes in the physical masses, leading to huge variations in the relic density 
especially when coannihilation processes are important or when annihilation occurs near 
a resonance. We will therefore take special care to disentangle the effects from kinematics 
and couplings. Indeed, as we will see, in many cases a large part of the phase dependence 
can be explained by changes in the masses of the involved particles. However, in some 
cases disentangling the kinematic effects will also lead to an enhancement of the phase 
dependence. 

The scenarios which we consider are the typical scenarios for neutralino annihilation: 
the mixed bino-Higgsino LSP that annihilates into gauge bosons, the rapid annihilation 
through a Higgs resonance, coannihilation with third generation sfermions, and finally a 
scenario with a mixed bino-wino LSP. The case of t-channel exchange of light sfermions 
is discussed together with the sfermion coannihilation. 

4.1 The mixed bino-Higgsino LSP 

We start with the case that all scalars except the light Higgs are heavy, Ms = mH+ = 
1 TeV. In this scenario we do not expect a dependence of the relic density on the phase 
of the slepton sector; we therefore set (pi = (pf In the real MSSM, a bino-like LSP with 
a mass of the order of 100 GeV needs a Higgsino admixture of roughly 25%-30% for its 
relic density to be within the WMAP range [14,66,67]. In terms of fundamental MSSM 
parameters this means Mi ^ fi. The main annihilation mechanisms then are XiX? ~^ WW 
and ZZ through t-channel chargino and neutralino exchange, as well as XiXi ~^ tt when 
kinematically allowed. The latter proceeds through s-channel Z or hi exchange. The LSP 
Higgsino fraction determines the size of the annihilation cross-section because it directly 
enters the XixfW^ and XiX^Z vertices. 

We perform a scan in the Mi-fi plane and display in Fig. la the region where the relic 
density is in agreement with the 2a WMAP bound, Eq. (1). In the real MSSM, when 
all phases are zero, only the narrow blue (dark grey) band is allowed. This band shifts 
slightly for negative values of /i (0^ = 180°). The onset of the ti annihilation channel 
appears as a kink. When allowing all phases to vary arbitrarily, keeping only those points 
for which all constraints are satisfied for at least one combination of phases, the allowed 
band becomes much wider, see the green (light grey) band in Fig. la. For a given Mi, the 
allowed range for fi increases roughly from Sfi ~ 10 GeV to 6n ~ 50 GeV. Since the eEDM 
constraint results in 0^ close to zero or 180°, this is actually mostly due to 0i. In fact the 
left boundary of the green band roughly corresponds to the contour of Qh"^ = 0.0945 for 
01 = 0^ = 180°. 

In Fig. lb, we show the relative mass difference between the lightest chargino and 
the LSP, Am-o-+/m^o = (m-± — m^o) I m^o , in the WMAP-allowed bands of Fig. la. As 
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Figure 1: (a) The 2a WMAP bands in the Mi-/i plane for tan/3 = 10, vnii+ = Ms = 
At = 1 TeV, for all phases zero (blue/dark grey band), for 0^ = 180° (or /i < 0) and 
all other phases zero (dashed red lines) and for arbitrary phases (green/light grey band), 
(b) The corresponding relative mass difference A = Am-o^+/m^o as function of m^o in 
the 2a WMAP band for all phases zero (blue/dark grey band) and for arbitrary phases 
(green/light grey band). 



a general feature, Am-o,-+/m^Q decreases with increasing m^o because the cross-sections 
for the pair-annihilation channels decrease, and coannihilation channels are needed in 
addition to maintain Qh"^ ~ 0.1. In the CPV-MSSM, however, for a given m^o much 
smaller mass differences can be in agreement with the WMAP bound than in the CP- 
conserving case. This is because, as we will discuss in more details below, the XiXi'W^^ 
couplings decrease with increasing 0i, so that additional contributions of coannihilation 
channels are required to maintain compatibility with WMAP. The phase dependence of 
the XiXi'W^''' couphngs is shown in Fig. 2 for Mi = 140 GeV and /i = 200 GeV. 

4.1.1 Below the tt threshold 

Figure 3 shows the WMAP band in the Mi-0i plane (0^ = 0^ = 0; = 0) for /i = 200 GeV 
and the other parameters as above. Also shown are contours of constant LSP mass, 
contours of constant LSP Higgsino fraction f^, as well as contours of constant cross- 
sections for the main annihilation channels. We can make several observations. First, 
the mass of the LSP increases with 0i. On the one hand this induces a decrease in 
the LSP pair-annihilation cross-sections. On the other hand, since the chargino mass 
is independent of 0i, the NLSP-LSP mass splitting is reduced, making coannihilation 
processes with xl (and also X2) more important. Second, the Higgsino fraction decreases 
with increasing 0i. The left- and right-handed XiXi W^"*" couplings feature the same phase 
dependence, see Fig. 2. The modification of the LSP couplings to gauge bosons leads to 
a decrease in the dominant XiXi ~^ WW/ZZ cross-sections and thus a higher value for 
the relic density. Here note that the phase dependence of fn (and hence of the couplings) 
is much more pronounced than that of the LSP mass. As a result, in Fig. 3a there is an 
almost perfect match between the 2a WMAP band and the band of 24% < fn "^ 28%. 
The small deviation at 0i ~ 40°-90° comes from the subdominant annihilation into Zhi 
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Figure 2: The X^Xi^~ couplings, 0^1 , Eq. (9), as function of 0i, for Mi = 140 GeV, 
/i = 200 GeV, (pf^ = (pt = 0, and the other parameters as in Fig. 1. The Higgsino fractions 
of the LSP, |A^i3p and jA^up, are also displayed. 

and hihi. The larger deviation at (pi ~ 150°-180° conies from coannihilations, c.f. Fig. 3b. 

These features also explain why in Fig. la nonzero phases only extend the WMAP- 
allowed range into the region where Qh"^ < 0.094 in the real case, and not into the one 
where Qh"^ is too large. Moreover, note that the phase dependence is large when annihi- 
lation into gauge bosons dominates, but weakened by contributions from coannihilation 
processes. 

To isolate the effect that comes solely from modifications in couplings, we display in 
Fig. 4 the variation of Qh"^ as function of 0i for constant m^o = 140 GeV. For comparison, 
the variation of ilh'^ for Mi = 164 GeV (leading to m^o ~ 140 GeV at (pi = 0) is also 
shown. In this example, flh"^ varies by a factor of 3 for constant Mi, and by a factor of 7 
for constant LSP mass. Overall, we find a phase dependence in Qh'^ of almost an order 
of magnitude for constant LSP mass. We emphasize that in this case of a mixed bino- 
Higgsino LSP, the dependence of masses and couplings on (pi work against each other, so 
that taking out the kinematic effects actually enhances the variation of Qh^. 

4.1.2 Above the ti threshold 

In the parameter range of Figs. 3 and 4 one is always below the tt threshold. We therefore 
consider in the next example a higher value of /i, such that XiXi ~^ tt contributes in the 
WMAP- allowed region. Analogous to Fig. 3, Fig. 5 shows the WMAP band in the Mi-0i 
plane for /i = 350 GeV together with contours of constant LSP mass, Higgsino fraction, 
and main (co) annihilation cross-sections. As can be seen in Fig. 5b, XiXi ~^ WW/ZZ and 
XiXi ~^ tt are of comparable importance, each contributing about 40%-50% to the total 
annihilation cross-section. The top-pair channel, proceeding through s-channel Z or hi 
exchange, shows a milder phase dependence (s-channel /i2,3 and t-channel stop exchange 
are negligible for this choice of parameters). Were it not for coannihilations, the lines 
of constant Qh'^ would again follow the lines of constant fn in Fig. 5a. However, since 
we are now dealing with a much heavier LSP, we also need a larger Higgsino fraction to 
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Figure 3: The 2a WMAP band (green/dark grey) in the Mi-(f)i plane for // = 200 GeV, 
tan/5 = 10, 171^+ = 1 TeV, At = 1.2 TeV, (j)t = (f)^ = 0. In the yellow (light grey) region 
Qh"^ is below the WMAP bound, and in the white region it is too large. Superimposed are 

125 and 140 GeV (dashed lines) and contours 
24% and 28% (dash-dotted lines). In (b) contours 



in (a) contours of constant LSP mass, m^o 



of constant LSP Higgsino fraction, fn 

Eq. (^ 
and the dashed lines for XiXi" coannihilation channels. 
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Figure 4: VlK^ as a function of 0i for the parameters of Fig. 3; the dashed line is for fixed 
Ml = 164 GeV, while for the full line Mi is adjusted such that m^o = 140 GeV. The 
green (grey) band shows the 2a WMAP-allowed range. 
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Figure 5: The 2a WMAP band (green/dark grey) in the Mi-0i plane for /i = 350 GeV, 
tan/3 = 10, mH+ = 1 TeV, At = 1.2 TeV, 0t = 0^ = 0. In the yellow (light grey) region, 
Qh"^ < 0.0945. Superimposed are in (a) a contour of constant LSP mass m^o = 292 GeV 
(dashed line) and contours of constant LSP Higgsino fraction fn = 29% and 37% (dash- 



dotted hues), and in (b) contours of constant {oJij.^kih'^) ^5 full lines: XiXi ~^ WW/ZZ, 



dashed lines: coannihilation with Xi and xi^ dotted lines 
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Figure 6: r2/i^ as a function of (^\ for the parameters of Fig. 5; the dashed line is for fixed 
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Figure 7: The 2a WMAP bands in the Mi-/i plane for mH+ = 500 GeV, At = 1.2 TeV 
and tan /? = 5 for all phases zero (blue/black bands) and for arbitrary phases (green/light 
grey bands). The other parameters are as in Fig. 1. In the yellow region, Qh"^ < 0.094. 
The dashed lines show contours of fi/i^ = 0.094, 0.128 when 0^ = 180°, 0i = 0. 

obtain the right relic density. At 0i = 0, this means fn — 29%-37%. This means in 
turn a smaller difference between Mi and n and hence a smaller Xi^Xi^ mass difference 
in the WMAP-allowed band as compared to the previous case. Therefore XiXt ^^id X1X2 
coannihilations are relatively more important. Since their cross-sections show the opposite 
phase dependence as those of the pair-annihilations, and since they are mainly determined 
by the mass difference, the overall phase dependence of the WMAP-allowed band is much 
weakened. 

In Fig. 6 we show the variation of Qh"^ as function of 0i for constant Mi = 314 GeV 
and for constant m^o = 292 GeV. We see that, taking out kinematic effects, for a relatively 
heavy bino-Higgsino LSP the phase dependence of Qh"^ is about a factor of 2. 

4.1.3 Lowering m.jj+ to 500 GeV 

Still keeping the sfermions heavy, we next lower the charged Higgs mass to mH+ = 
500 GeV. This leads to rapid annihilation via s-channel Higgs exchange when m^o ^ 
mh.-/2 ^ 250 GeV. In this case only a very small coupling of the LSP to one of the 
Higgses (/i2, /is) is necessary, so /i can be large and the LSP still annihilates efficiently 
even if it is dominantly bino. 

In Fig. 7 we display the WMAP-allowed regions in the Mi-/i plane for both the real 
MSSM and the CPV-MSSM. As before, in the CPV-MSSM, the allowed region corre- 
sponds to the points in the Mi-fi plane for which there exists at least one choice of 01, 
0^, (pt, 4>i for which all constraints are satisfied. One clearly sees the two very narrow 
bands of the so-called Higgs funnel at m^o ?» 171^+/ 2 and large fi. The impact of the 
phases on the relic density in the funnel region will be discussed in the next section. 
Outside this region, i.e. for small /i, we observe as before a significant widening of the 
band consistent with WMAP for nonzero phases. Furthermore, for 0^ 7^ 0, lower values 
of /i give ?n^± > 103.5 GeV, consistent with the LEP constraint [68] on charginos. This 
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Figure 8: Cross section times relative velocity, crt;(xiXi ~^ W~^W ), as function of 0i for 
p = 10 GeV, Ml = 220 GeV, /x = 260 GeV, At = 1 TeV and the other parameters as in 
Fig. 7. The long-dashed line is the xf exchange, the short-dashed line the hi exchange, 
and the full line the total contribution. Here p is the center-of-mass momentum and 

V = 2p/m^o. 

widening of the WMAP band is again due to shifts in couplings and masses as we have 
discussed above. Again, the widening is into the region where Vth'^ < 0.0945 in the CP- 
conserving case. Here note that for Mi = 200 GeV the extra allowed region to the right 
of the blue band corresponds to 0^ = 180°. 

The shifts in masses are specially relevant in the region around the Higgs funnel. For 
example at Mi = 270 GeV, the band of allowed values for /i increases from 15 to 160 GeV 
when allowing for arbitrary phases. Here one is still close enough to the resonances to 
have dominant annihilation through Higgs exchange, and small changes in the /i2,3 masses 
have a large effect on Vth"^. Furthermore, because the couplings of the LSP to the heavy 
Higgses can be suppressed, in the CPV-MSSM the LSP can have a much larger Higgsino 
component as compared to the real case and still be in agreement with WMAP. 

A priori one could think that in this region, where several channels contribute to the 
LSP annihilation, the impact of the shifts in couplings could be amplified by interference 
effects, leading to a significant impact on the effective annihilation cross-section. Although 
we do find interference effects, they have little influence on the WMAP-allowed bands 
shown in Fig. 7. In fact, contrary to what was originally reported in [69], the interference 
effects between the s-channel Higgs and t-channel chargino exchange diagrams for XiX? ~^ 
W^W~ are destructive, so that an enhancement of the cross section in one channel is 
largely cancelled by the other channel. An example for such an interference between the 
t-channel Xi ^^^ s-channel hi exchange diagrams is shown in Fig. 8. 

4.2 Annihilation through Higgs 

In the Higgs sector, nonzero phases, in particular 0j, can induce scalar-pseudoscalar mix- 
ing as well as important changes in the masses. One can therefore expect large differences 
between the real and complex MSSM in the Higgs-funnel region. 

At vanishing relative velocity, t" — > 0, neutralino annnihilation through s-channel scalar 



ifi 



exchange is p-wave suppressed; the annihilation proceeds strictly through pseudoscalar ex- 
change. Nevertheless when performing the thermal averaging, the scalar exchange cannot 
be neglected altogether. In the MSSM with real parameters it can amount to (9(10%) 
of the total contribution. In the presence of phases, all the neutral Higgs bosons can 
acquire a pseudoscalar component (that is gj^ -o -o 7^ 0) and hence significantly contribute 
to neutralino annihilation even at small v. There is a kind of sum rule that relates the 
couplings squared of the Higgses to neutralinos. Therefore, for the two heavy eigenstates 
which are in general close in mass, we do not expect a large effect on the resulting relic 
density from Higgs mixing alone. A noteworthy exception occurs when, for kinematical 
reason, only one of the resonances is accessible to neutralino annihilation. That is for 
example the case when mh2 < 2m^o ~ rra/13. 

For the analysis of the Higgs funnel, we choose 

Ml = 150 GeV, tan/5 = 5, Mg = 500 GeV, At = 1200 GeV. (27) 

We consider the cases of small and large Higgsino mass parameter, /i = 500 GeV and 
fi = 1-2 TeV, leading to small and large mixing in the Higgs sector respectively for 
(pt 7^ 0. In both cases the LSP is dominantly bino. As mentioned above, allowing for 
nonzero phases not only affects the neutralino and Higgs couplings but also their physical 
masses. Since the relic density is very sensitive to the mass difference Arriy^Oj^. = mhi—2m^a 
[13,23], it is important to disentangle the phase effects in kinematics and in couplings. 
As we will see, a large part of the huge phase effects reported in Ref. [49] can actually be 
attributed to a change in Am^o^.. 

4.2.1 Small Higgs mixing 

We fix yU = 500 GeV so that there is small Higgs mixing. Moreover, we set 0^ = to 
avoid the eEDM constraint, and discuss the dependence on 0i and (f)f 

For real parameters and mH+ = 340 GeV, we have m^o = 147 GeV, rrih^ = 331.5 GeV, 
m/t3 = 332.3 GeV, and Qh"^ = 0.11. In this case, /i2 is the pseudoscalar. The LSP annihila- 
tion channels are then characterized by the /i2 branching fractions, giving predominantly 
fermion pairs, bb (78%) and rf (10%), with a small contribution of Zhi (7%). When we 
vary the phases of At and/or Mi , we observe large shifts in the relic density. 

First consider varying the phase (pt, which affects the Higgs masses and mixings 
through loop effects. In this scenario with relatively small //, the scalar-pseudoscalar 
mixing never exceeds 8%. In Fig. 9a we plot the band that is allowed by WMAP in the 
TT^H+^'Pt plane. One can see that the lower and upper WMAP bounds correspond to the 
contours of Am^o^^ = 36.2 and 38.6 GeV respectively, with only 4% deviation. So the 
main effect of (pt can be explained by shifts in the physical masses and position of the 
resonance. 

We next vary the phase 0i, keeping (pt = 0. This changes the neutralino masses and 
mixing, and hence also the neutralino-Higgs couplings, Eq. (18). For m//+ = 340 GeV, 
when increasing 0i, the relic density drops, see Fig. 9b. This is because the mass of 
the neutralino increases slowly, resulting in a smaller Am^o^^- Adjusting m^o or mh^ 
(by changing Mi or mH+) such that the mass difference stays constant, we find rather 
that the relic density increases with 0i. This can be readily understood from the phase 
dependence of the couplings of /i2,3 to the LSP, shown in Fig. 10. For 01 = 0, the coupling 
of /i2 is predominantly pseudoscalar and h^ almost purely scalar while for (pi = 90°, it is 
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Figure 9: The 2a WMAP bands (green/dark grey) in the mH+-(j)t and mH+-(j)i planes for 
the parameters of Eq. (27) and /i = 500 GeV; 0i = in (a) and 0^ = in (b). Contours 
of constant mass differences Am^o^^ = ''^/i2 ~ ^.m^o are also displayed. In the yellow (light 
grey) region, Qh"^ is below the WMAP range. 



hs that has a large pseudoscalar coupling. Therefore for (pi = 0, /i2 exchange dominates 
with a large cross section while for 0i = 90° one gets about equal contributions from /i2 
and hs, although with a smaller overall cross section. When increasing 0i further (up 
to 180°), /i2 exchange again dominates, however with a coupling to neutralinos smaller 



0. Thus one needs a smaller mass splitting Am^Of^^ for ilh to 



by 30% than for 0i 

fall within the WMAP range, see Fig. 9b. Moreover, for large 0i there is also a sizeable 
contribution from XiXi ~^ hihi with a constructive interference between s-channel h^ 
and t-channel neutralino exchange. In Fig. 11 we show the variation of Vth"^ with 0i while 
keeping Atti^o^^ fixed. The maximum deviation which comes purely from modifications 
in the couplings can reach 70%. 



4.2.2 Large Higgs mixing 

We next discuss the case of large mixing in the Higgs sector, which is achieved for large 
values of /i = 1-2 TeV. For this purpose we concentrate on the (f)t dependence. It has to 
be noted that here, for large (pt ~ 90°, rather light /i^ and large /iv4j, the EDM constraint 
is not satisfied by simply setting (p^ = 0. One has to either allow for a small (p^ ~ 0{1°), 
or appeal to cancellations due to light sleptons with masses of few xlOO GeV. We choose 
the latter solution, imposing 0^ = 0. Adjusting the selectron parameters such that 
de < 2.2 X 10~^'^ecm has only a (9(1%) effect on the relic density in our examples. 

In Figure 12 we show the WMAP-allowed regions in the mfj+-Mi plane for this choice 
of parameters and maximal phase of At {cpt = 90°). The regions for which Qh"^ falls within 
the WMAP band are shown in green (dark grey), and those for which Qh^ is too low in 
yellow (light grey). In addition, the positions of the WMAP-allowed strips for (pt = are 
shown as dashed lines. In the CP-conserving case, h^ is a pure pseudoscalar and /i2 a 
pure scalar, while for (pt = 90° it is just the opposite and /i2 is dominantly pseudoscalar. 

For /i = 1 TeV, Fig. 12a, the mass splitting between /i2,3 is about 10 GeV for (pt = 90°, 
as compared to about 2 GeV for (pt = 0. Masses and pseudoscalar contents, H^^, of /i2,3 are 
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Figure 10: The scalar (full lines) and pseudoscalar (dash lines) neutralino Higgs couplings, 
g ^0^0, as a function of (pi for 171^+ = 340 GeV, and the other parameters as in Fig. 9b. 
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Figure 11: Qh"^ as a function of 0i for mH+ = 340 GeV and the value of Mi adjusted so that 
Am^o^2 stays constant (full line). The other parameters as in Fig. 9b {rrih^ = 332.5 GeV) 
For comparison, the variation of Qh"^ for fixed Mi = 150 GeV is shown as a dashed line. 
The green (grey) band corresponds to the 2a WMAP range. 
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Figure 12: The WMAP-allowed bands in the m//+-Mi plane for (a) /i = 1 TeV and 
(b) /i = 2 TeV with 0^ = 90°, 0i = and the other parameters given by Eq. (27). In the 
narrow green (dark grey) bands 0.0945 < Qh'^ < 0.1287, while in the yellow (light grey) 
regions Qh"^ < 0.0945. The positions of the WMAP bands for (pt = are shown as dashed 
lines. 



depicted in Fig. 13 as function of (pf Here note that it is /i2, i-e. the state which changes 
from scalar to pseudoscalar with increasing 0j, which shows the more pronounced change 
in mass; the crossovers of 50% scalar-pseudoscalar mixing of /i2,3 occur at (pt ~ 15° and 
145°. For Ml values up to 250 GeV, we therefore find in both the CP-conserving and the 
CP- violating case two narrow bands where 0.094 < Qh"^ < 0.129. For 0j = (and also for 
0j = 180°) both these bands are mainly due to pseudoscalar /13 exchange, with one band 
just below and the other one above the pseudoscalar resonance. For 0^ = 90° the situation 
is different: in the lower WMAP-allowed band the LSP annihilates through the scalar hs, 
with the pseudoscalar /i2 not accessible because tti/ij < 2m^o ~ mh-^^, while in the upper 
band both /i2 and h^, contribute (with /i2 exchange of course dominating). In between the 
two WMAP-allowed green bands one is too close to the pseudoscalar resonance and Qh'^ 
falls below the WMAP bound; this holds for both (pt = and (pt = 90°. The positions 
of the WMAP-allowed bands for 0j = and (pt = 90° are not very different from each 
other. Still the difference in the relic density between (pt = and (pt = 90° is typically 
a factor of a few in the WMAP-bands, and can reach orders of magnitudes at a pole. 
For Ml > 250 GeV and (pt = 90°, one enters the region of coannihilation with stops, 
leading to a vertical WMAP-allowed band. For (pt = 0, the ti is 55 GeV heavier, so the 
stop coannihilation occurs only at Mi ~ 305 GeV (for (pt = 180° on the other hand, 
rut^ ~ 230 GeV and coannihilation already sets in at Mi ~ 200 GeV). 

For /i = 2 TeV, Fig. 12b, there is an even stronger CP-mixing of /i2,3 and the mass 
splitting between the two states becomes ~ 45 GeV for (pt = 90°. The pseudoscalar 
contents are similar to those in Fig. 13 (r.h.s. plot) with the 50% cross-over at (pt ~ 20°. 
Moreover, because the LSP has less Higgsino admixture, one has to be closer to resonance 
to obtain the right relic density. As a result, the scalar and pseudoscalar funnels become 
separated by a region where Qh"^ is too large [50]. In fact both the /i2 and h^ exchange each 
lead to two WMAP-allowed bands, one above and one below the respective resonance. 
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Figure 13: Masses and pseudoscalar content of h^ and /i3 as function of 0( for mii+ = 
335 GeV, yU = 1 TeV, and the other parameters given by Eq. (27). The hght Higgs has a 
mass of rrih^ ^ 117 GeV and a pseudoscalar content of < lO"'^. 



250 




70 80 
Mi (GeV) 



100 



Figure 14: Same as Fig. 12a but for hght neutrahnos and Higgs bosons. The vertical red 
line indicates the LEP bound on the xl mass, the horizontal dashed line that on rrih^. 



For the h^ (scalar) exchange, however, these two regions are so close to each other that 
they appear as one line in Fig. 12b. This is in sharp contrast to the CP-conserving case, 
0i = 0, where the scalar and pseudoscalar states are close in mass, hence leading to only 
two WMAP-allowed bands. These are again shown as dashed lines in Fig. 12b and origin 
dominantly from the pseudoscalar resonance, the scalar resonance being 'hidden' within. 
We also display in Fig. 14 the WMAP-allowed bands for /i = 1 TeV focusing on 
the region of small neutralino and Higgs masses. Here we see clearly three specific Higgs 
annihilation funnels, each delimited by two narrow bands where Qh"^ is within the WMAP 
bound. The band corresponding to hi exchange is also found in the CP-conserving case. 
However, in the CP-conserving case the LEP limit [35] on the Higgs mass, rriho > 114 GeV, 
requires mfj+ > 210 GeV, while in the CP-violating case the limit is much lower, about 
rrihi > 85 GeV [70] at tau/S = 5 and nit = 175 GeV, corresponding to mH+ > 130 GeV. 
For mj|/+ < 190 GeV, the bands corresponding to /i2 and /13 exchange are clearly separated 
because here the mass splitting between the two states is large, about 20-36 GeV. This 
is to be contrasted with the real case, where the H/A mass splitting is much smaller and 
annihilation through the pseudoscalar always dominates. 
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Figure 15: The 2a WMAP bands (green/dark grey) in the 'mH+-(j)t plane for yU = 1 TeV 
and mH+ ~ 330-340 GeV in (a) and mH+ ~ 305 GeV in (b). Contours of constant mass 
differences Am^Oj^^ are also displayed. In the yellow (light grey) region, Qh'^ is below the 
WMAP range. 






Let us now examine the dependence on 0j in more details. For this we fix /i = 
1 TeV. For vanishing phases, agreement with WMAP is found for mH+ — 332-334 GeV. 
This value is lower than in the scenario with small Higgs mixing because the Higgsino 
fraction of the LSP is smaller, so one needs to be closer to the Higgs resonance. For 
(pt y^ we have a large scalar-pseudoscalar mixing and hence a stronger dependence of 
Qh"^ on (pt. For (pt = 0, h^ is the pseudoscalar and gives the dominant contribution 
to neutralino annihilation while for (pt = 90° /i2 is the pseudoscalar, hence giving the 
dominant contribution. Consequently in Fig. 15, agreement with WMAP is reached for 
25 GeV with hi = /13 at 0* = and 180°, and hi = /i2 at 0* = 90°. 

When twice the LSP mass is very near the heaviest Higgs resonance, one finds another 
region where the relic density falls within the WMAP range. This is shown in Fig. 15b 
(corresponding to the phase dependence of the lower WMAP-allowed band in Fig. 12a). In 
the real case one needs 171^+ = 305 GeV, giving a mass difference Am^Ofj^^ = —1.5 GeV. 
Note that annihilation is efficient enough even though one catches only the tail of the 
pseudoscalar resonance. For the same charged Higgs mass, the mass of h^ increases when 
one increases (pt, so that neutralino annihilation becomes more efficient despite the fact 
that /13 becomes scalar-like and g?o-Of^ decreases. When 0t ~ 75°-90°, the coupling 
g^o-Of^ becomes very small and one needs Am^Of^^ = 0-1.5 GeV to achieve agreement 
with WMAP, see Fig. 15b. Here we are in the special case where rrih^ < 2m^o ~ "Ti/ig, so 
that only h^ contributes significantly to the relic density. 

We can isolate the phase dependence of Qh'^ due to the scalar-pseudoscalar mixing 
by keeping the distance from the /13 pole constant. For constant values of Am^Of^^ = 
— 1.5 GeV we get an increase in Qh^ relative to the 0j = case by almost an order of 
magnitude. This is however far less than the huge shifts of several orders of magnitude 
found for fixed values of mH+ when a Higgs pole is passed. 
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4.3 Bino-like LSP and light sfermions 

In the CP-conserving MSSM, light sfermions can significantly contribute to reducing the 
relic density to a value which is in agreement with WMAP, in particular in the case of a 
bino-like LSP. The relevant processes are XiXi ~^ ff via t-channel sfermion exchange, as 
well as coannihilation with sfermions that are close in mass to the LSP. In the CPV-MSSM 
with large phases, the sfermions of the first and second generation need to be heavy to 
avoid the EDM constraints. The third generation is however much less constrained and 
can be light. We therefore consider in this section the cases of light staus and light stops. 

4.3.1 Light staus 

We choose a scenario where the LSP is mostly bino and fix tan/9 = 10, /i = 600 GeV, 
mH+ = Ms = Af^t = 1 TeV. Moreover, we take M^^ = 220 GeV and M^^ = 240 GeV, 
so that staus are relatively light, nif-^^ = 212.1 GeV, and can contribute to the neutralino 
annihilation. Note that for this choice of parameters there is a large mixing in the stau 
sector, driven by yU.tan/3. The eEDM constraint is avoided by setting (pfj, = (pt = 4>i = 0. To 
obtain a relic density in agreement with WMAP one has to rely on stau (co) annihilation 
channels. For this aim, the mass difference Amj^o^^ = rrif-^ — m^o must be small, in our 
example about 4-6 GeV. For Mi = 210 GeV and vanishing phases, we obtain m^o = 
208 GeV, Am^Of^^ =4.1 GeV and Qh"^ = 0.102. The main channels in this case are 
Xin -^ IT (32%), x?ri ^ Zt (10%), nn -^ tt (26%) and nf* ^ 77 (12%). 

It is well known that the mass difference is the key parameter in case of coannihilations. 
We therefore expect large shifts in Qh"^ for nonzero phases, resulting from small changes 
in the masses. Figure 16a shows the WMAP-allowed region in the Mi-(j)i plane for the 
scenario given above. As can be seen, the WMAP band matches almost perfectly with 
the contours of constant mass difference, Am^Of_^ = 3.7 and 5.6 GeV. When we adjust 
the parameters of the stau sector to keep a constant mass difference while varying 0i, the 
relic density stays constant within few percent, 6Q/Q < 5%. An analogous behaviour is 
found for nonzero (pi. 

For very light neutralinos and sleptons, annihilation into lepton pairs through t-channel 
slepton exchange can be efficient enough to achieve Qh"^ ~ 0.1. In the mSUGRA model, 
this is often called the 'bulk' region. Owing to the LEP limit of mj > 90-100 GeV [68] 
(depending on the slepton flavour and chirality/mixing and on mj—m^o), and because the 
t-channel slepton contribution scales as mlo/'m^f, the 'bulk' is squeezed into a small region 
of slepton masses of about 100 GeV. To investigate this case in the CPV-MSSM, we lower 
the stau parameters to M^^ = 135 GeV, M^^ = 150 GeV and A^ = 100 GeV, (pi = 0. 
This gives rrif, = 107.8 GeV, ruf^ = 182.1 GeV and m^^ = 135.7 GeV. The WMAP band 
in the Mi-(p plane for this scenario is shown in Fig. 16b. For m^o < 100 GeV, that is 
up to Ml ~ 100 GeV, XiXi ~^ t~^t~ is the dominant process. For heavier neutralinos, 
X?"?! coannihilation dominates. Agreement with WMAP is achieved for larger Xi^^i mass 
differences than in the previous example. Since coannihilations are less important, we flnd 
a stronger phase dependence which is not completely determined by Am^Of^, see Fig. 16b. 
When keeping the masses constant, the maximal variation in Qh'^ due to 0i is about 15%. 
Last but not least note that s-channel Z and Higgs exchange is negligible in this scenario; 
t-channel exchange of f2 however does contribute and there is in fact a strong destructive 
interference between the fi and 'f2 exchange diagrams. 
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Figure 16: The 2a WMAP bands (green/dark grey) in the Mi-cpi plane for m^^ ~ 212 GeV 
in (a) and rrif^ ^ 108 GeV in (b). The other parameters are given in the text. In the 
yellow (light grey) region, Qh'^ is below the WMAP bound. Contours of constant mass 



difference Am^o^ are shown as dashed lines. 
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Let us lower the stau masses parameters even further, close to the experimental limit 

130 GeV, Mf = 140 GeV, we get m^, = 98.2 GeV, 



0. 



For M^3 



while keeping 0, 

771^2 = 175.9 GeV and nij^^ = 124.6 GeV. The WMAP-allowed regions in the Mi-0i plane 
for this case are shown in Fig. 17. The almost horizontal bands of annihilation through 
the light Higgs {rrih^ = 118 GeV) for very light Xi; as well as the stau coannihilation 
region at Mi > 90 GeV are clearly visible. The peculiar feature is that for large phase, 
01 > 110°, there appears a new region, connecting the light Higgs funnel and the stau 



~o~o 



coannihilation strip, where Qh is within the WMAP bound. In this region, xix 



lAl 



TT 



completely dominates, and the phase dependence originates from the Xi^i,2T couplings. 
Again, there is an important interference between the fi and f2 exchange diagrams. 

That this new 'stau bulk' region is generic can be seen in Fig. 18. Here we plot the 
WMAP-allowed bands in the Mq-Mi plane, with Mq = M^^ = 0.9Mi^. For vanishing 
phases, the light Higgs funnel and f coannihilation regions appear as narrow disconnected 
strips. For arbitrary phases 0i, 0^, 0f, 0/, these strips are much wider; especially the light 
Higgs funnel becomes a band instead of a narrow strip. Furthermore, the Higgs funnel 
is connected to the f coannihilation region by the t-channel stau exchange region, which 
appears as a horizontal band at Mq ~ 130 GeV. 



4.3.2 Light stops 

To discuss the case of a light stop, we fix Mq^ 
800 GeV, n = Ms = mH+ = 1 TeV, and tan/5 = 
the eEDM constraint. Moreover, we choose At = 



Mu. 



= 500 GeV 
5. We again fix 
1 TeV and (hf = 



Mn, 



= 450 GeV, 
bf^ = to easily avoid 
180° to obtain a light 



ti, m^^ = 243.5 GeV. Since /i is large we are again in a scenario with a bino LSP. A relic 



density in agreement with WMAP is found for Mi 
real case. The main channels are XiXi ~^ tt (21% 



215 GeV 



' Xi 



X% -^ thi (57%), and x?i~i 



214 GeV) in the 

gt 



(17%). Note that Xi^^i coannihilation dominates although the mass difference is much 
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Figure 17: The 2a WMAP bands (green/dark grey) in the Mi-0i plane for mf^ ~ 98 GeV. 
The other parameters are given in the text. In the yellow (light grey) region, VLh"^ is below 
the WMAP bound. Also shown are contours of constant mass differences Atti^o^^ = 

m/ij — 2m^o (dash-dotted lines). 
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Figure 18: The 2a WMAP bands in the Mq-Mi plane, with Mq 

blue (black) bands are for vanishing phases, the green (medium grey) one for arbirary 

phases, 0i, (pt = 4>i- The pink (dark grey) region is excluded by the requirements m^± > 
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LSP. 
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Figure 19: The 2o" WMAP bands (green/dark grey) for the stop coannihilation scenario, 
(a) in the Mi-0i plane for rrif:^ = 243.5 GeV and (b) in the At-4't plane for m^o = 
210.8 GeV. The other parameters are given in the text. In the yellow (light grey) regions, 
Qh'^ < 0.0945. The medium grey region in (b) is theoretically excluded. Also shown are 
contours of constant mass difference Am^o^^ (blue dashed). The full black line in (a) is a 
contour of constant (o"f)(xi^i —^ thi). 



larger than in the case of Xi^i coannihilation. 

The phases that can play a role here are 0i and (pf As we have already observed 
in other coannihilation scenarios, the relic density is extremely sensitive to the mass 

However, in the stop coannihilation scenario we also 



TTlf 



Xi 



difference Am^or 

observe some important effects due to the phase dependence of the couplings. 

First we vary only (pi and show in Fig. 19a the WMAP-allowed band in the Mi- 
(pi plane. We find the WMAP band does not match the contours of constant mass 



difference Am^o^^ 



mz 



Xi 



A much larger mass difference is required at 



180° 



(Am^Oi^ = 34.2-38.6 GeV) then at (^i = (Am^oj^ = 28.8-31.0 GeV). The reason for this 
is an increase of both the XiXi ~^ tt as well as the Xi^i ~^ thi cross-sections with the phase 
01. For the coannihilation process, the phase dependence is enhanced by a constructive 
interference between the t-channel ti and the s-channel top exchange diagrams. 

To investigate also the dependence on (pt, we now fix Mi = 212 GeV and 0i = 
(m^o = 210.8 GeV) and plot in Fig. 19b the WMAP-allowed band in the At-(pt plane. 
The other parameters are as above. As before, agreement with WMAP is found only 
for a narrow band in which x^ti — > thi dominates. Although in this band the ti mass 
is constant within 10 GeV, there is no exact match between the contours of constant 
Qh"^ and Am^o^^. Rather, at (pt = 0, 90° and 180°, agreement with WMAP requires 
36 GeV, 32 GeV and 30 GeV, respectively (each about ±2 GeV). 

When keeping the mass difference constant, we can single out the phase dependence 
of Qh!^ that comes solely from changes in the couplings. This is shown in Fig. 20, where 
we plot Qh"^ as a function of (pi for (pt = 180° (dashed line) and as a function of (pt for 
01 = (dash-dotted line), each time keeping Am^ot-^ = 32 GeV constant by adjusting 



Am^oi^ 



either Mi or At. As can be seen, in either case Qh^ varies between 0.08 and 0.14. 
comparison, the variation of nh'^ with 0i for Mi = 212 GeV {At = 1 TeV, (pt = 180°) 
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Figure 20: Qh"^ in the stop coannihilation scenario as a function of 0i (dashed hne) and 
as a function of 0t (dash-dotted hne) for constant mass difference Am^oi^ = 32 GeV. For 
comparison, the variation of Qh"^ with (pi for fixed Mi = 212 GeV is also shown. See text 
for details. 

varying mass difference, is also shown. 

4.4 Relaxing the gaugino GUT relation: the bino/wino scenario 

We now consider relaxing the universality condition amongst the first two gaugino mass 
parameters and treat Mi and M2 as two independent parameters. As we want to examine 
specifically the role of the wino component we fix /^ = 1 TeV and mH+ = 1 TeV, then the 
Higgsino component of the LSP is small and the annihilation near Higgs resonance not 
possible. Because we choose sfermions to be heavy, the annihilation into fermion pairs is 
suppressed. The LSP can hence only pair-annihilate into gauge bosons. 

Choosing M2/M1 < 2 increases the wino component of the LSP. Since the neutralino- 
chargino-W^ coupling has also a term proportional to Ni2, one could think that the 
bino-wino scenario is quite similar to the mixed bino-Higgsino case discussed in Sec- 
tion 4.L However, when the parameters are set such that fw = |A'i2p becomes sizable, 
say around 10%, the mass difference Am-o^+ becomes small (few GeV), and the coanni- 
hilation channels are so important that the relic density is well below the WMAP limit, 
unless Mi,2 ~ C(l) TeV. 

In fact, in the mass range which is interesting for collider searches, the LSP has to be 
still overwhelmingly bino, and the relic density is completely dominated by coannihilation 
channels involving x?xi^, X?X25 X2Xi^, X2X2> xfxf-^ Final states involve a variety of 
channels from gauge boson pairs to fermion pairs of all three generations. Even when 
letting all phases vary, the relic density falls within the WMAP range only for a very 
narrow range of parameters; for a given value of Mi the viable range of M2 varies by only 
2-3 GeV. For tan/9 = 10, for instance, the allowed band can roughly be parameterized as 



M2 ~ (2.03 X lO-^Mi^ + I.O2M1 + 9) ± 2 



(2^ 



with Ml given in GeV. This is typical for scenarios that are dominated by coannihilation 
■^This scenario was considered in Rcf. [71] for negative values of Mi, that is (f>i = 180°. 
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Figure 21: The 2cr WMAP-allowed band (green/dark grey) in the Mi-(j)i plane for M2 = 
200 GeV, tan/3 = 10, fj, = Ms = mn+l TeV, (f>t = 4>fj. = 0- Superimposed are contours 
of constant x1 masses (blue dashed lines). In the yellow (light grey) region, Qh'^ is below 
the WMAP range. 



in the sense that only a narrow range of mass difference between the NLSP and the 
LSP is allowed. Furthermore, considering the large number of contributing channels, 
the phase dependence in each individual channel tends to be "softened". For example 



10, we have Am 



Xi Xi 



21.2 GeV and 



for M2 = 200 GeV, Mi = 179 GeV, tan/3 
Qh"^ = 0.121; the dominant channels are xtxt ~~* f f altogether contributing around 40% 
of the total effective annihilation cross-section. In Fig. 21, we display the WMAP allowed 
range in the Mi-{/)i plane for </)^ = 0. The contours of constant Vth^ basically follow the 
contours of constant mass difference Ar?7,-o^+, as expected when coannihilation channels 
are dominant. Only for (pi > 90° there is a small increase in Qh'^ due to shifts in couplings, 
introducing a small gap between the contours of constant mass difference and those of 
constant Qh"^. For a given Am-o^+, the maximal deviation from the case of vanishing 
phases reaches AQ/Q ^ 25%. 



5 Summary and Conclusions 

We have performed the first complete computation of the relic density of neutralino dark 
matter in the CPV-MSSM, including in a consistent manner phases in all annihilation 
and coannihilation processes. Moreover, we have presented a comprehensive study of the 
typical scenarios that predict a relic density in agreement with WMAP. Since CP phases 
do not only change the sparticle and Higgs mixings but also the masses, we have taken 
care to disentangle effects from kinematics and couplings. 

A priori one could think that taking out the effects which come from changes in the 
masses would diminish the phase dependence of Qh^. For processes, for which the mass 
difference is the most important parameter, i.e. for coannihilations or annihilation near a 
pole, this is indeed the case. One the other hand, we have found several examples where 
the phase dependence of masses and couplings work against each other, and taking out 
the kinematic effects actually enhances the variation in Qh"^. This happens, for instance. 
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in the case of a mixed bino-Higgsino LSP, where we have found effects of almost an order 
of magnitude from modifications in the couphngs due to nonzero CP phases. In the case 
of annihilation through s-channel Higgs exchange, for small scalar-pseudoscalar mixing, 
we have found cases where for fixed MSSM parameters Qh^ goes down with increasing 
01, but when keeping the mass difference between the pseudoscalar-like Higgs pole and 
2m^o constant, Qh"^ actually goes up. This effect is of the order of 50% and dominantly 
due to changes in the scalar-/pseudoscalar-type coupling of the LSP. Much larger effects 
have been found for large scalar-pseudoscalar mixing. In the peculiar case that only 
one resonance is accessible to neutralino annihilation, we have found order-of-magnitude 
variations in Qh^ due to changes in the pseudoscalar content of this resonance. Moreover, 
in some cases we have found large interference effects in the most important annihilation 
or coannihilation channels. Such interferences do, however, not necessarily lead to a large 
variation in Qh^. 

We have also considered scenarios with a bino-like LSP and light sfermions (stops and 
staus), and studied the phase dependence in the regions where a) annihilation into // 
dominates, b) coannihilations dominate, as well as the intermediate region where c) both 
t-channel exchange and coannihilation are important. In the stop-coannihilation region, 
when x5^i —>■ thi is the dominant process, the variation in Qh"^ can be about due to 
changes in the Xi^i^i couplings. There is moreover a constructive interference between 
the t-channel ti and the s-channel top exchange diagrams. Another peculiar feature arises 
for (very) light staus: here we have found a region where for large phase 0i the annihilation 
into TT alone can be efficient enough to obtain agreement with the WMAP bound. This 
region does not occur for zero phases. 

We emphasize that even in scenarios which feature a modest phase dependence once 
the kinematic effects are singled out, the variations in Qh"^ are comparable to (and often 
much larger than) the ~ 10% range in fi/i^ of the WMAP bound. Therefore, when 
aiming at a precise prediction of the neutralino relic density from collider measurements, 
it is clear that one does not only need precise sparticle spectroscopy — one also has to 
precisely measure the relevant couplings. As we have shown in this paper, this programme 
certainly includes the determination of possible CP phases. For zero phases, in a bino 
scenario with light sleptons, one may be able to infer Qh"^ of the LSP with roughly the 
WMAP precision (~ 15%) from LHC measurements [72,73]. At the ILC, one expects 
to achieve much higher precisions, allowing for a prediction of AQ/Q of the level of few 
percent in case of a bino LSP annihilating through light sleptons [74] and in the stau 
coannihilation scenario [75,76]. In the bino-Higgsino scenario, AQ/Q ~ 15% may be 
achieved [74]. Whether similar precisions can be reached in the CPV-MSSM requires 
careful investigation. 

To emphasize the need to determine as completely as possible the underlying parame- 
ters of the model, we stress again that we have found examples where in the CPV-MSSM 
the relic density of the LSP could be quite different as compared to that in the MSSM 
with vanishing phases. Simply from a precise measurement of part of the mass spectrum 
one could be lead to wrongly conclude that, for instance, the model does not give a good 
dark matter candidate, or else that some significant thermal production is necessary to 
explain the observed number density. At the same time it is important to note that a 
computation of (av) at next-to-leading order will be necessary to achieve the required 
precision in the prediction of Qh"^, see [77-79]. Last but not least we stress that in addi- 
tion direct or indirect detection of the CDM candidate will be indispensable to pin down 
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the dark matter in the Universe. 

Finally we remark that we have not taken into account the constraint arising from 
b -^ s'-f. The supersymmetric corrections to the branching ratio for b ^ s'~f depend 
mostly on the squark and gaugino/Higgsino sector as well as on the charged Higgs. In the 
real MSSM, large corrections are found at large values of tan/5, which we do not consider 
in this paper. Large corrections might also arise in scenarios with a light charged Higgs. 
A detailed study of the impact of this measurement in the CPV-MSSM is left for a future 
work. 

Note added: While this paper was in preparation, the WMAP collaboration published 
new results corresponding to three years of data taking. The new WMAP+SDSS com- 
bined value for the relic density of dark matter is QcDuh"^ = O.lll^o^oii ^^ -'-'^ l^^]- This 
is only slightly below the value used in this paper, so our conclusions do not change with 
the new data. 



Acknowledgements 

We thank T. Gajdosik and W. Porod for comparisons in the CPV-MSSM and J.S. Lee 
for his help with CPsuperH. We also thank S.Y. Choi for discussions. This work was 
supported in part by GDRI-ACPP of CNRS and by grants from the Russian Federal 
Agency for Science, NS-1685.2003.2 and RFBR-04-02-17448. S.K. is supported by an 
APART (Austrian Programme for Advanced Research and Technology) grant of the Aus- 
trian Academy of Sciences. A. P. acknowledges the hospitality of CERN and LAPTH 
where some of the work contained here was performed. 

A Interaction Lagrangian 

We here give the relevant sparticle interactions with other sparticles and SM particles in 
the CPV-MSSM. The Higgs boson interactions with sparticles and particles can be found 
in [53]. 



A.l Neutralino— neutralino— Z 

Neutralinos couple to Z bosons via their Higgsino components: 
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One can also write Eq. (29) as 
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A. 2 Chargino— chargino— Z, 7 

The interaction of two charginos with electroweak gauge bosons is 

l^x^x^z = j^Z;^^^ (o;^ P^ + 0:f P^ X- - e A, ^Yxt (32) 

with i,j = 1, 2, and 

O'i = -V,,V;,-lV,,V*, + 6,^smHw, (33) 

0[f = -U:,U^, -lU:,Uj, + 6,^ sin' Ow. (34) 

A. 3 Neutralino— chargino— W 

The neutrahno-chargino-W interaction is described by {i = 1, 2; j = 1, ...,4) 

^xxw = 9 W-^^ r {O}. Pl + Of, Pj,) xt + 9 W+^ ^ {Of* P^ + Of* P^) x] (35) 

with 

Of = N^,V:,~j=^N^,V:, and Of = N*^U,, + ^^N*,U,,. (36) 

A. 4 Neutralino— fermion—sfermion 

The sfermion interaction with neutrahnos is [i = 1,2; j = 1, ...,4) 

%x° = 9fif[,PR + h{^PL)x',fL + 9fihi^PB + flt,PL)xU~R + ^-^- 

= 9fUPR + hl^P^y^f, + h.c. (37) 



where 



a 



f _ 

ij JLj-^m I '"Rj 



/i4 + 4,4, (38) 

b!, = hi^Ri + fi^Ri. (39) 

The couphngs f[ ^ and /i{^ ^ are 

/i- = -72 (^^-2 + i tan ^H'iV,!) , 4. = ^ (iV^.2 - I tan O^Nj,) , (40) 

/4. = ^ tan^H/iV;, , fl = -^ tan^H^iV;i , (41) 

4 = -^* ^.-4 = /^S , <■ = -K N,3 = h% (42) 
for stops and sbottoms, and 

4. = ^(tane^% + iV,2), (43) 

4. = -V2tan^H/iV;i, (44) 

/^li, = -K Nfs = hi) (45) 
for staus. In more general terms, 

fl = -V2 ((e; - 4) tan^^iV,! + liN,2) , 4 = v^e^ tan^n^iV;, . (46) 
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A. 5 Chargino— fermion— sfermion 

The sfermion interaction with charginos is {i,j = 1,2) 

^/'/x± = 9ui-UjiPj, + KV;,PL)xtdL + 9uiKUj2PR)xtdR + 

= gu{lt^Pj, + ^Pj}xtd, + gd{ir,Pn + ^PL)xpu, + \i.c. (47) 

where u {u) stands for up-type (s)quark and (s)neutrinos, and d (d) stands for down-type 
(s) quark and charged (s)leptons. The couphngs / and k are 

4 = -VjiRl + K Vj2RL , 4 = -u,,Ri + K u,A , (48) 



-* — /)L n* R* k^ 

for stops and sbottoms and 



k\^ = h U*Rl, , 4 = ht V*Rli , (49) 



l^ = -Vj, , /[,. = -Uj.Rl, + K Uj2Rl^ , (50) 

k^ = KU*,, kj^ = (51) 

for staus and sneutrinos. 

A. 6 Sfermions with gauge bosons 

The sfermion interaction with photons is the same as in the CP-conserving case: 

Cjj^ = -tee J A, if I d^ f\ + /* d^ /^) = -tee J 5,, A, f* d^ f, . (52) 

The interaction with Z bosons is given by 

^ffz = '-^Z.iCUld^h + CnfRd^h) 

■' ■' cos Uw 



7W 



cos 9 



w 



(c,r(m; + CrrLrI;) z, i; o^ h m 



with Cj;^ j^ = lli^ R~ ^f sin^ 6w- Note that there is only a phase dependence for i j^ j. 
The interaction with W bosons is given by 



V2 



Cjj,^ = -^ {w;^iid^h + w; h\ dn^) 



= -^ (^L4* w;^ i* d^ k + R^^P^ w;: h] d^ i,) (54) 

taking t bW as an example for simplicity. The corresponding FejTiman rules are obtained 
from 

Ad''B = A{d,B)-{d,A)B ^ f*d^f^ = i{k^ + k^Y (55) 

where k^ and k- are the four-momenta of /^ and /• in direction of the charge flow. 
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